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Abstract 

Using fragmentation functions for charged pions, charged kaons, and (anti)pro- 
tons recently extracted from experimental data of e^e~ annihilation at the Z-boson 
resonance and at centre-of-mass energy ^/s = 29 GeV, we perform a global study of 
inclusive charged-hadron production in pp, jp, and 77 collisions at next-to-leading 
order in the parton model of quantum chromo dynamics. Comparisons of our re- 
sults with pp data from CERN SppS and the Fermilab Tevatron, 7p data from 
DESY HERA, and 77 data from CERN LEP2 allow us to test the universality of 
the fragmentation functions predicted by the factorization theorem. Furthermore, 
we perform comparisons with e'^e~-annihilation data from LEP2 so as to test the 
scaling violations predicted by the Altarelli-Parisi evolution equations. 

PACS numbers: 13.60.Le, 13.65. -hi, 13.85.Ni, 13.87.Fh 



1 Introduction 



In the framework of the QCD-improved parton model, the inclusive production of single 
hadrons is described by means of fragmentation functions (FFs), D^{x,n'^). The value 
of D^{x,jj?) corresponds to the probability for the parton a produced at short distance 
to form a jet that includes the hadron h carrying the fraction x of the longitudinal 
momentum of a. Unfortunately, it is not yet possible to derive the FFs from first principles, 
in particular for hadrons with masses smaller than or comparable to the asymptotic 
scale parameter, A. However, given their x dependence at some energy scale /i, the 
evolution with /i may be computed perturbatively in QCD using the timelike Altarelli- 
Parisi equations [|I|]. This allows us to test QCD quantitatively within one experiment 
observing single hadrons at different values of centre-of-mass (CM) energy -y/i (in the case 
of e^e~ annihilation) or transverse momentum pt (in the case of scattering). Moreover, 
the factorization theorem guarantees that the D^(x,/i^) functions are independent of the 
process in which they have been determined, and represent a universal property of h. 
This enables us to make quantitative predictions for other types of experiments as well. 

During the last five years, the experiments at the CERN Large Electron-Positron 
Collider (LEPl) and the SLAC Linear Collider (SLC) have provided us with a wealth 
of high-precision information on how partons fragment into low-mass charged hadrons. 
The data partly comes as light-, c-, and 6-quark-enriched samples without or 
with identified final-state hadrons (vr^, K^, and p/p) or as gluon-tagged three-jet 

samples without hadron identification P,pi][. This new situation motivated us to update, 
refine, and extend a former analysis performed by two of us together with Binnewies [p!T| , p!^ 
by generating new leading-order (LO) and next-to-leading-order (NLO) sets of vr^, i^^, 
and p/p FFs fl^.Q By also including in our fits ti^, K^, and p/p data (without flavour 



separation) from SLAC PEP [Tj], with CM energy = 29 GeV, we obtained a handle 
on the scaling violations in the FFs, which enabled us to determine the strong-coupling 
constant. We found af\Mz) = 0.1181 ± 0.0085 at LO and a^^^Mz) = 0.1170 ± 0.0073 
at NLO |]T3|. These results are in perfect agreement with what the Particle Data Group 
currently quotes as the world average, a^f\Mz) = 0.1181 ± 0.002 ||16[. Our strategy was 



to only include in our fits LEPl and SLC data with both flavour separation and hadron 
identification p,^,||J^ gluon-tagged three-jet samples with a fixed gluon-jet energy 0, and 
the vr^, K^, and p/p data sets from the pre-LEPl/SLC era with the highest statistics 



and the finest binning in x |]T^. Other data served us for cross checks. In particular, 
we probed the scaling violations in the FFs through comparisons with vr^, K^, and p/p 
data from DESY DORIS and PETRA, with CM energies between 5.4 and 34 GeV fT^. 
Furthermore, we tested the gluon FF, which enters the unpolarized cross section only at 
NLO, by comparing our predictions for the longitudinal cross section, where it already 



FORTRAN subroutine which returns the values of the D'^(x,ii'^) functions for given values of x 
and may be downloaded from the URL http://www.desy.de/~poetter/kkp.html or obtained upon 
request from the authors. 

^Contrary to what is stated in the respective figure captions, the squares pertaining to the upmost 
curves shown in Figs. 2-4 of Ref. @ refer to OPAL data §]. 
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enters at LO, with available data PJT8|. Finally, we directly compared our gluon FF with 
the one recently measured by DELPHI in three-jet production with gluon identification 
as a function of x at various energy scales fi |10]. All these comparisons led to rather 
encouraging results. We also verified that our FFs satisfy reasonably well the momentum 
sum rules, which we did not impose as constraints on our fits. 

The purpose of this paper is twofold. On the one hand, we wish to extend our previous 
tests of scaling violations [|l^] to higher energy scales by confronting new data on inclusive 
charged-hadron production in e~^e~ annihilation from LEP2 [22,23], with y/s ranging from 
133 GeV up to 189 GeV, with NLO predictions based on our FFs. On the other hand, 
we wish to quantitatively check the universality of our FFs by making comparisons with 
essentially all available high-statistics data on inclusive charged-hadron production in 
colliding-beam experiments. This includes pp data from the UAl p^p5|,p6| and UA2 



27,28| Collaborations at the CERN Super Proton-Antiprotron Synchrotron (SppS) and 
from the CDF Collaboration E^] at the Fermilab Tevatron, 'jp data from the HI 



and 

ZEUS Collaborations at DESY HERA, and 77 data from the OPAL Collaboration 
at LEP2. 

Very recently, alternative sets of NLO FFs for vr^, K"^ 



T9[] , and charged hadrons [ p^ , pO[ 



have become available. They are based on different collections of experimental data and on 
additional theoretical assumptions. In Ref. |T9|, power laws were assumed to implement 
a hierarchy among the valence- and sea-quark FFs. In Ref. pO| , the renormalization and 
factorization scales were identified and adjusted according to the principle of minimal 
sensitivity . In order to estimate the present systematic uncertainties in the FFs, it is 
useful to compare NLO predictions for inclusive charged-hadron production consistently 
evaluated with the three new-generation FF sets [|l^,|l^,^. 

This work is organized as follows. In Section |^, we present comparisons of our NLO 
predictions for inclusive charged-hadron production with e^e 
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pp 



data from LEP2 

data from SppS |l||,|6|,|3>il and the Tevatron pf, 7p data from HERA iO|,|ll, and 77 
data from LEP2 . In Section |^, we compare our NLO FFs with those of Refs. p!9|j20| by 
studying inclusive charged-hadron production in e+e~ annihilation at various CM energies 
for the case of quark-fiavour separation. Our conclusions are summarized in Section ^ In 
the Appendix, we list convenient parameterizations of the x and n"^ dependences of our 
7r±, K^, and p/p FFs at LO and NLO pi. 



2 Global Analysis of Collider Data 

Throughout this paper, we work at NLO in the MS renormalization and factorization 
scheme with n/ = 5 active quark fiavours. We choose the renormalization scale /i and 
the factorization scale Mf to be /i = Mf = ^a/s in the case of e^e~ annihilation and 
/i = Mf = C,Pt otherwise, where the dimensionless parameter ^ is introduced to determine 
the theoretical uncertainty due to scale variations. As usual, we allow for variations of ^ 
between 1/2 and 2 around the default value 1. As for the parton density functions (PDFs) 
of the proton, we employ set CTEQ5M provided by the CTEQ Collaboration |]53|. As 
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for the photon PDFs, we use the set by Aurenche, Fontannaz, and Guillet (AFG) |]3^ as 
default and those by Gliick, Reya, and Vogt (GRV) and by Gordon and Storrow (GS) 
for comparison. As for the charged-hadron FFs, we take the sum of our tt^, i^'^, and 



pp FFs [|l3l as default. In Section y, we also use the charged-hadron sets by Kretzer (K) 
||T9| and by Bourhis, Fontannaz, Guillet, and Werlen (BFGW) EO] for comparison. We 



evaluate ah'^\ij) from the two- loop formula [jT6l with the value of A^ that belongs to the 



respective FF set, namely A^ = 213 MeV for Ref. |13|, A- 



Ag^ = 215 MeV for Ref. |2g. These values are in the same ball park as the one pertaining 



MS 



MS 



168 MeV for Ref. 19 , and 



MS 



to the CTEQ5M || PDFs, A^ = 226 MeV. The considered photon PDFs [0,11,11 



are all implemented with A^ 



- 200 MeV, which corresponds to A^ = 137 MeV for a 
4.5 GeV. The NLO cross sections for the partonic subprocesses 



6-quark mass of nn, 

e~^e^ —^c + X,ab^c + X,'yb^c + X, and 77 c + X, where a, b, c denote 
partons (quarks and gluons) and X represents the residual final-state partons, were found 

1^,^, respectively. The NLO formalism for calculating 
annihilation, hadroproduction in pp collisions, and 



in Refs. 0, !§, Hgg, and 
the hadronic cross sections for e 



photopro duct ion in e^p and e"'"e~ collisions is described in Refs. |12,|2,|3,H|, respectively. 
In the following, we always consider the sum of positively and negatively charged hadrons. 
While the experimental data are sampled in bins, we present our theoretical predictions 
as continuous distributions. This should be justified as long as the bins are narrow or the 
data is corrected for finite bin widths. 



2.1 



Annihilation 



In Ref. [|T^ , we already tested the scaling violations in our FFs, predicted by the timelike 
/i^ evolution equations, by making comparisons with e~^e~ data of inclusive vr^, K^, and 
p/p production in the energy range from ^/s = 5.2 GeV to 44 GeV [l^. The advent of 
LEP2 data on e+e" ^ + X with = 133, 161, 172, 183, and 189 GeV, taken by 

p3| , enables us to gradually extend these tests towards high 



DELPHI 22 and OPAL 



energies. This data is given as the normalized differential cross section [1 / a tot) da / dx as a 
function of the scaled momentum x = 2p/ -y/i, where p is the absolute value of the three- 
momentum of the final-state hadron. In Fig. |I], it is compared with our NLO predictions 
for ^ = 1/2, 1, and 2. 

We observe that the scale variation of the latter is rather small, about ±1% at x = 
0.1 and ±10% at x = 0.8, although compensating effects on the renormalization scale 
dependence only set in at next-to-next-to-leading order. In order to assess the agreement 
between experimental data and theoretical predictions quantitatively, we evaluate the 
values per degree of freedom, x^p, for ^ = 1/2, 1, and 2. As in Ref. |jl3|, we select the 
data points with x > 0.1. The resulting Xdf values and the respective numbers of data 
points are summarized in Table 0. The Xbf values are of order unity, just as for the data 
used in our fits 
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22,23 



We thus conclude that inclusion of the LEP2 data 
would have had little impact on the resulting FFs. For y/s 
are best described by the central scale choice ^ = 1, while the data at higher energies 



in our fits 
133 and 161 GeV, the data 
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tend to prefer ^ = 1/2. From the nice overall agreement we conclude that the predicted 
scaling violations are confirmed by experiment. 



with 1 < \y\ < 1.8 |28 
= 630 GeV 



2e 



2.2 Hadroproduction in pp Collisions 

The differential cross section Ed^a/d^p = {l/2TTpT)d'^a/dy dpx of pp + X was 

measured at SppS and the Tevatron as a function of pt at various values of ^/s in rapidity 
(y) intervals constrained by the respective detector geometry, mostly in the central region. 
Specifically, UAl provided data at = 200 [H], 500 540 H, 630 HU, and 900 GeV 
II, UA2 at = 540 GeV p3,|8[, and CDF at = 630 and 1800 GeV p9[. This data 
is compared with our NLO predictions for ^ = 1/2, 1, and 2 in Figs, ^j-^ The respective 
Xqp values are presented in Table |I]. They are evaluated using the data points with 
Pt > /io/^5 where /xq = \/2 GeV is the starting scale of our FFs [|^, so that backward 
evolution is avoided. For ^/s values of order 500 GeV, the scale variation ranges from 
about ±60% at Pt = 3 GeV to about ±25% at pr = 25 GeV, and at a given value of 
Pt it decreases with increasing values of ^/s. The data is generally well described by our 
central predictions, both in normalization and shape, even at low values of pt, where 
corrections from beyond NLO and nonperturbative effects are expected to be significant. 
However, the UAl data sets with = 200 |5[ and 630 GeV |2§ and the UA2 data set 

^ tend to prefer the scale choice ^ = 1/2. For the UAl data at 

this was also noticed in Ref. [^. In that paper, it was also found 
that the overall agreement between this data and the NLO prediction for ^ = 1/2 does 
not improve if the normalization factor for the gluon FF is modified within the tolerance 
allowed by the original fit and the other FF parameters are in turn adjusted so as to 
optimize the fit. As for the UA2 data with 1 < \y\ < 1.8 [^] and the UAl data at 
^/s = 630 GeV |j26|, we should also remark that the experimental situation is somewhat 
ambiguous, since the UAl data at ^/s = 540 GeV |2^ and the CDF data at ^/s = 630 GeV 
clearly favour ^ = 1. 

At first sight, one might be worried by the relatively large Xdf values obtained for 
^ = 2. However, we have to bear in mind that, according to the prescription underlying 
Table |l|, they are evaluated from all data points with pt > /io/2 ~ 0.71 GeV. On the 
one hand, the low-pr data from the hadron colliders come with very small errors. On the 
other hand, perturbation theory is known to break down in the \ow-pT limit. Therefore, 
the Xdf values for ^ = 2 in Table |I] should be considered as exploratory. In order to 
have a more conservative basis to judge the agreement between experiment and theory 
and to facilitate the comparison between the cases ^ = 1/2, 1, and 2, we present in 
Table || the corresponding Xdf values obtained by selecting the data points with pt > 
2/io ~ 2.83 GeV, independent of ^. This happens at the expense of reducing the number 
of data points contributing to the evaluation of the Xdf values. In particular, the UA2 
data with \y\ < 0.63 ||2^ is excluded altogether, and just two out of the 22 CDF data 
points at y/s = 630 GeV are left. Comparing Table ^with Table ^ we observe that 



the Xdf values for ^ = 1 and 2 are considerably reduced, especially for the UAl data at 



v/s = 630 GeV ^ and the CDF data at ^/s = 1800 GeV 
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2.3 Photoproduction in e^p Collisions 



HI (ZEUS measured the differential cross section dPa/dy dp\ of inclusive charged- 
hadron photoproduction in collisions of 27.5 GeV positrons (26.7 GeV electrons) and 
820 GeV protons, with ^ = 300 GeV (296 GeV), as a function of pt for |?/| < 1 
(—1.2 < y < 1.4) in the laboratory frame. Here, y is taken to be positive in the proton 
beam direction. In this experiment, the scattered lepton was tagged, so that the virtuality 
Q'^ of the Weizsacker- Williams photon was limited to be < 0.01 GeV^ (0.02 GeV^), 
and the 7p invariant mass W was taken to be in the interval 165 < W < 251 GeV 
(167 < W < 194 GeV). HI also measured the y distributions da/dy for px > 2 and 
3 GeV. In Figs. this data is compared with our NLO predictions for ^ = 1/2, 1, and 
2. The scale variation approximately amounts to ±35% at = 3 GeV and to ±15% at 
Pt = 12 GeV. The corresponding Xbf values are given in Tables |l] and @. As in the case 
of the pp data, they are evaluated using the data points with pt > fio/C, and pr > ^fiQ, 
respectively. In the case of the y distribution for pt > 2 GeV, the NLO prediction for 
^ = 1/2 also involves energy scales fi < fio and should, therefore, be taken with a grain of 
salt. The corresponding Xdf value in Table |l| is marked by an asterisk. We observe that 
the Pt and y distributions from HI are nicely described by the theoretical predictions 
for ^ = 1, as for both normalization and shape, with Xdf ^ 1- the other hand, the 
ZEUS data |^ prefer ^ = 2 or even larger values of ^. 

In Fig. ^ we also investigate the dependence of the NLO prediction on the photon 
PDFs. We observe that the y distributions evaluated with the GRV set have a slightly 



larger normalization than those evaluated with the AFG |33] set, but they essentially have 
the same shape. By contrast, the y distributions based on the GS |3^ set are appreciably 
steeper. However, we should keep in mind that the these differences are smaller than the 
theoretical uncertainty due to scale variations. 



2.4 Photoproduction in e+e Collisions 

OPAL 1^ measured the differential cross section of e^e^ e~^e~h^ +X via two-photon 
collisions at LEP2, with average CM energy i/i = 166.5 GeV, in the 77-invariant-mass 
interval 10 < W < 125 GeV and in the three subintervals 10 < < 30 GeV, 30 < W < 
55 GeV, and 55 < 1^ < 125 GeV. The scattered electrons and positrons were antitagged, 
the maximum scattering angles being 6' = 33 mrad. Specifically, they presented the pt 
distributions da/dpT for \y\ < 1.5 and the \y\ distributions da/d\y\ for pt > 1.5 GeV. 
Notice that the y distributions are symmetric w.r.t. y = 0, so that da/d\y\ = 2da/dy. 
This data is compared with our NLO predictions for ^ = 1/2, 1, and 2 in Figs. ^ and p!0| , 
respectively. In the case of the calculation for 10 < ly < 125 GeV, the scale variation 
approximately amounts to ±20% at px = 3 GeV and to ±5% at px = 15 GeV. The 
corresponding Xbf values are given in Tables ^ and 0. As in the case of the pp and 7p 
data, they are evaluated using the data points with px > /io/^ and pr > 2/io, respectively. 
In Fig. |lO[ the NLO predictions for ^ = 1/2 involve energy scales /i < /iq, as is indicated 
in Table |l| by asterisks, but they are nevertheless shown for illustration. From the steep 
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fall-off of the pt distributions in Fig. p| it follows that the y distributions in Fig. [10 



are 



greatly dominated by the low-pr range, where the theoretical uncertainty due to scale 
variations is considerably larger than the experimental errors and nonperturbative effects 
may be significant. Thus, the pt cut-off in Fig. |10| is likely to be too small for perturbation 
theory to be meaningful. We observe that the px and y distributions are nicely described 
by the NLO predictions for ,^ = 1, as for both normalization and shape. 

In Figs. |ll] and |12|, we illustrate how our default NLO predictions (evaluated with 
the AFG photon PDFs for ^ = 1) for 10 < < 125 GeV shown in Figs. | and 0, 
respectively, are built up by their direct, single- resolved, and double- resolved components. 
As expected, the double-resolved channel dominates at low values of pt, while the direct 
channel wins out at large values of px- However, one should keep in mind that this 
decomposition depends on the scale choice and the factorization scheme. In Figs. |Tl] and 
|12|, we also analyze the theoretical uncertainty related to the photon PDFs. Similarly to 
Fig. ^ the NLO predictions based on the AFG and GRV sets only differ in their 



overall normalizations. By contrast, the pt and y distributions evaluated with the GS [36 



set exhibit a somewhat different shape. There is actually a minor caveat attached to the 
latter: the interval of pt integration involves values below the starting scale /iq = GeV 



of the GS ||36[ set. The comments concerning the validity of perturbation theory made 



in connection with Fig. |TO|also apply to Fig. |T2|. We conclude that, except in the low-pr 
range, the theoretical uncertainty due to the photon PDFs is insignificant, of order 10% 
for Pt > 4: GeV. 



3 Comparisons with Other FF Sets 



In order to compare the three NLO FF sets that were released earlier this year |[13|jl9|j20| , 
we consider the differential cross section {l/atot)dcr/dx of inclusive charged-hadron pro- 
duction from the fragmentation of light (uds), c, b, and all quarks in e~^e~ annihilation at 
a/s = 29, 91.2, and 189 GeV. We compare the resulting NLO predictions for ^ = 1 with 
PEP ra, LEPl mm, SLC m, and LEP2 [Ell data. In contrast to Ref. [BOl, we refrain 



from scale optimization. Specifically, we consider the all-flavour samples from DELPHI 
§, SLD §, and OPAL |2|], the light-quark samples from TPC |5| and DELPHI |,|^, 
the c-quark samples from TPC ||4^ and DELPHI [^, and the 6-quark samples from TPC 
45| and DELPHI 



Notice that the DELPHI analysis of Ref. Q is based on data 
collected during the years 1991-1993, while the one of Ref. [0 is based on 1994 data. The 
outcome of these comparisons is presented in Figs. [l3Hl5|, which refer to ^/s = 29, 91.2, 
and 189 GeV, respectively, and Table 0. 

In order to properly assess these comparisons, it is useful to recall which of these data 
samples actually entered the fits of Refs. |T^Jl^,^ . The flavour-specific data samples from 
TPC were used in Ref. 
1994 data from DELPHI 



T|], the 1991-1993 data from DELPHI § in Ref. gg], the 



in Refs. ig,|2D[f| and the SLD § data in Refs. 1T|JT|,|2D 
each of the three analyses [Il3| , p!9| , p0 



For 

the data samples which were not used as input are 



^The light-quark sample from Ref. was not used in Ref. ||2C 
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marked in Table ^ by asterisks. We did not include the flavour-specific data samples from 
TPC in our fit [|l^ because they are unpublished. Furthermore, we only incorporated 
the most recent data from DELPHI 0. In contrast to Refs. [12,19,20|, we did not make 
use of the ALEPH data, for reasons explained in Ref . [T^ .[ 
From Figs. [U 



T5|, we observe that the three NLO FF sets 



T|,|T|,gg yield very simi- 
lar all-flavour and light-quark results, except for the small-x region {x<,0.1), where the 
QCD-improved parton model is known to be unreliable, and the large-x region {x ^0.8), 
where the experimental data is less constraining anyway. On the other hand, there are 
appreciable differences in the c- and 6-quark contributions also at intermediate values of 

13 



X. 



In the c-quark case, the K 19 



131 and BFGW 



ones, 



result undershoots the KKP 
while the KKP result undershoots the K [|l3] and BFGW |^ ones in the 6-quark 
case. The latter feature may be understood by observing that the 6-quark FFs of Ref. [|1^ 
were determined so as to fit the more recent DELPHI [0 data, whereas, for example, the 
BFGW I^D]] analysis used the older DELPHI data as input. As may be seen from 



Fig. |TJ, the 6-quark data samples of these two measurements 0,0 are not compatible 
with each other. This was also observed in Ref. [^. Thus, it is impossible to simultane- 
ously describe both 6-quark data samples with Xdf values of order unity, as is also 
evident from Table Leaving aside this 6-quark anomaly, we conclude from Table ^ that 
the three NLO FF sets |lT3|JT9| , pO| are compatible with each other and nicely describe the 
considered data samples H,[7|J^,^ , with Xdf values of order unity, even in those cases 
where the latter were not included in the fits. 

In this section, we compared the cross sections of charged-hadron production from the 
fragmentation of light, c, b, and all quarks resulting from the three NLO FF sets 
rather than the individual FFs at some particular value of /i^. Comparisons of the latter 
kind were performed in Refs. [|r^,0 taking into account the old BKK sets rather 
than the new KKP [|1^ ones, and appreciable differences were found, in particular at large 
values of x. We feel that, for the time being, such comparisons are of limited usefulness, 
especially as far as the individual u-, d-, and s-quark FFs are concerned. These FFs are 
not separately constrained by available e~^e~ data, and different theoretical assumptions 
were made in Refs. IIT3|,n!9|,|20 



4 Conclusions 

Although the FFs yU^) characterizing the hadronization of partons a into light 

hadrons h are genuinely nonperturbative objects, they possess two important proper- 
ties that follow from perturbative considerations within the QCD-improved parton model 
and are amenable to experimental tests. Firstly, their dependence on the factorization 
scale /i is supposed to be determined by the timehke Altarelli-Parisi equations [|1|. Sec- 
ondly, the factorization theorem predicts that they are universal in the sense that they 

"'Due to a misinterpretation of Refs. H,^, the respective Xbf values for CT''(all), (T^{uds), a'''{c), and 
(t''(6) quoted in Table 1 of Ref. ||l3| are incorrect. The correct numbers read 48.1, 11.1, 7.61, and 104, 
respectively. 
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only depend on the fragmenting parton a and the inclusively produced hadron h, but 
not on the details of the process from which they were determined. In this paper, we 
tested these two properties at NLO for the tt^, K^, and p/p FFs recently fitted [|1^ to 
experimental data of e'^e~ annihilation from LEPl [^^0, SLC [§, and PEP On the 
one hand, we confronted new data on inclusive charged-hadron production from LEP2 
2^,|23|, with CM energies ranging from 133 GeV to 189 GeV, with our NLO predictions. 



On the other hand, we performed a global NLO analysis of essentially all high-statistics 
data on inclusive charged-hadron production in colliding-beam experiments, including pp 
scattering at SppS [^^,^,^,^ and the Tevatron |^ , 7p scattering at HERA |30|j31 
and 77 scattering at LEP2 



32 



In all cases, we found reasonable agreement between 
the experimental data and our NLO predictions as for both normalization and shape, 
as may be seen from Figs. |l]-0 and Tables ^ and 0. The majority of the data sets are 
best described with the central scale choice ^ = 1. Exceptions include the UAl data sets 
200 |2§ and 630 GeV |2§ and the UA2 data set with 1 < 
1/2, as well as the ZEUS data, which favours ^ = 2. 



with a/s 
which prefer ^ 



\y\ < u 



and the UA2 data set with 

However, if we 

estimate the theoretical uncertainty due to unknown corrections beyond NLO by varying 
^ between 1/2 and 2, as is usually done, then it is justified to state that all the considered 
data sets agree with our NLO predictions within their errors. We hence conclude that 
our global analysis of inclusive charged-hadron production provides evidence that both 
the predicted scaling violations and the universality of the FFs are realized in nature. 

As is well known, the gluon FF enters the prediction for the unpolarized cross section 
da/dx of inclusive hadron production in e~^e~ annihilation only at NLO, while at LO it 
only contributes indirectly via the /i^ evolution. In order to nevertheless have a handle 
on it, we included in our fits |jl3| experimental data on gluon-tagged three-jet events from 
LEPl [P]. Furthermore, we checked that our predictions for the longitudinal cross section. 



where it already enters at LO, agree well with available data P,|18|. On the other hand 



the gluon FF is known to play a crucial role for pp, •yp, and 77 scattering at low values 
of pt- Thus, the comparisons performed here provide another nontrivial test of the gluon 
FF. 

As we have seen in Figs. ^|-^ |^, and the theoretical uncertainty of the NLO predic- 
tions due to scale variations significantly decreases as pt increases. In order to perform 
more meaningful comparisons, it would, therefore, be desirable if pp, 7p, and 77 experi- 
ments extended their measurements out to larger values of px- Furthermore, in order to 
render such comparisons more specific, it would be useful if these experiments provided 
us with separate data samples of n'^, K^, and p/p hadrons. 

We also estimated the current systematic uncertainties in the FFs by comparing our 
NLO predictions for flavour-tagged inclusive charged-hadron production with those ob- 
tained from two other FF sets which appeared during this year |lT9|pO| . Apart from a 
difference in the fe-quark FFs at medium to large values of x, which may be traced to the 
incompatibility of two underlying 6-quark-specific data samples all three FF sets 

P^ , P!P| , PD[ mutually agree within the present experimental errors. We conclude that the 
determination of NLO FFs from global fits to experimental data has reached a level of 
precision which is comparable to the one familiar from similar analyses for PDFs. 
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A Parameterizations 

For the reader's convenience, we present here simple parameterizations of the x and /i^ 
dependences of our LO and NLO sets of vr^, K^, and p/p FFs [jl3[. As in Ref. [0, we 
introduce the scaling variable 

'-'(^s/aS)- ^ 

In LO (NLO), we have A^^ = 88 MeV (213 MeV) We use three different values for 
yUo, namely H 

{^2 GeV, if a = u,d,s,g 

m(r/e) = 2.9788 GeV, if a = c . (2) 

m(T) = 9.46037 GeV, if a = b 

This leads to three different definitions of s. For definiteness, we use the symbol Sc for 
charm and Sf, for bottom along with s for the residual partons. We parameterize our FFs 
as 

D(x,/i2) = Ax"(l-x)^(^l + ^) (3) 

and express the coefficients N, a, /3, and 7 as polynomials in s, Sc, and Sb- For s = Sc = 
Sb = 0, the parameterizations agree with Eq. (2) of Ref. |]13| in combination with the 



appropriate entries in Table 2 of that paper. The charm and bottom parameterizations 
must be put to zero by hand for Sc < and Sb < 0, respectively. 

We list below the parameters to be inserted in Eq. (^. The resulting parameterizations 
describe the properly evolved FFs to an accuracy of 10% or less for fj,Q < fi < 200 GeV 
and 0.1 < X < 0.8. Deviations in excess of 10% may occur outside these regions. 

LO FFs for (7r+ + vr-): 

Z^f (x,/i2) = Df (x,/i2): 



N = 0.54610 - 0.22946s - 0.22594r + 0.21119r 
a = -1.46616 - 0.45404s - 0.126845^ + 0.27646s^ 
P = 1.01864 + 0.95367s - 1.09835s^ + 0.74657s^ 
7 = -0.01877s + 0.02949s^ (4) 
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N = 22.2815 - 20.8125s - 11.5725s^ + 15.5372s^ 
a = 0.12732 + 0.23075s - 2.71424^^ + 1.724568^ 
(3 = 6.13697 + 2.18849s - 5.04475s^ + 3.29117s^ 
7 = 0.09044s - 0.075895^ 

N 
a 

7 

N 
a 

7 

N 
a 

7 

NLO FFs for (7r+ + 7r-): 

N = 0.44809 - 0.13828s - 0.069515^ + 0.01354s^ 
a = -1.47598 - 0.30498s - 0.01863s^ - 0.12529s^ 
(3 = 0.91338 + 0.64145s + 0.07270s^ - 0.16989s^ 
7 = 0.07396s - 0.077575^ 

N 
a 
/3 
7 
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= 8.75500 - 9.32277sc + 1.80600s^ + 2.02179s^ 
= -0.38611 - 0.41190sc - 0.48496s^ + 0.42525s^ 
= 5.61846 + 0.74035sc - 0.64929s^ + 0.66788s^ 
= 0.06652sc - 0.05531s^ 



= 0.31147 - 0.19319sfe - 0.10487s^ + 0.18824s^ 
= -1.92993 - 0.44692sb - 0.08271s^ + 0.30441sj| 
= 3.47086 + 0.79775S6 - 0.28091s^ + 0.39504s^ 
= -0.04887sb + 0.03212s^ 



= 6.04510 - 6.61523s - 1.649785^ + 2.68223s^ 
= -0.71378 + 0.14705s - 1.084235^ - 0.431825^ 
= 2.92133 + 1.48429s + 1.32887s^ - 1.78696s^ 
= 0.23086s - 0.29182s^ 



= 16.5987 - 18.3856s + 2.44225r + 2.13225s^ 
= 0.13345 + 0.22712s - 0.83625s^ + 0.38526s^ 
= 5.89903 - 0.16911s + 0.59886s^ - 0.25630s^ 
= -0.18619s + 0.87362s^ 



N = 6.17173 - 4.82450sc - 1.30844s^ + 1.95527s^ 
a = -0.53618 - 0.27879sc - 0.513375^ + O.lOQOOs^ 
(3 = 5.60108 + 0.83571sc - 1.15141s^ + 0.77027s^ 
7 = 0.09268sc- 0.11267s^ 

N = 0.25944 - 0.11499S6 + 0.03733s^ - 0.18028s^ 
a = -1.98713 - 0.35858S6 + 0.22277s^ - 0.66413s^ 
P = 3.52857 + 0.72303S6 + 0.46260s^ - 0.99235s^ 
7 = -0.02701S6 - 0.02089s^ 

N = 3.73331 - 3.16946s - 0.47683s^ + 0.70270s^ 
a = -0.74159 - 0.51377s - 0.197055^ - 0.179175^ 
f3 = 2.33092 + 2.03394s - 0.507648^ - 0.08565s^ 
7 = 0.09466s- 0.10222s2 

LO FFsfor {K+ + K'): 

N = 0.25937 - 0.10502s + 0.005725^ - 0.002695^ 
a = -0.61925 + 0.09956s + 0.073895^ - 0.00070s^ 
p = 0.85946 + 0.57965s + 0.263978^ - 0.127648^ 
7 = 0.15303s + 0.14807^2 

N = 5.38115 - 3.05084s - 1.10056s^ + 1.31207s^ 
a = -0.00321 - 0.25889s - 0.184945^ + 0.13994s^ 
(3 = 3.07632 + 1.13745s - 0.90413s^ + 0.56581s^ 
7 = 0.05141s - 0.00697s^ 

N = 5.18266 - 3.48519sc - 1.00982s^ + 1.17996s^ 
a = -0.17751 + 0.02309sc - 0.61327s^ - 0.03532s^ 
P = 4.30306 + 1.00547sc - 0.51779s^ + 0.20683s^ 
7 = 0.13514sc-0.17778s^ 



12 



N = 1.57044 - 1.78340sb + O.STlOOs^ + 0.15469s^ 
a = -0.84143 - 0.43448s6 - 0.05314s^ - 0.36621s^ 
P = 6.01488 + 0.72953S6 - 0.64433s^ + 0.92351sg 
7 = 0.01024S6 - 0.06160s^ 



N = 0.02862 - 0.02113s + 0.00389s^ + 0.009015^ 
a = -2.94091 + 0.66881s - 0.29670s^ + 0.20574s^ 
(3 = 2.73474 - 0.58222s + 0.04329s^ + 0.78033s^ 
7 = 0.03586s - 0.012205^ 



NLO FFs for {K+ + R-): 



N = 0.17806 - 0.10988s - 0.02524s^ + 0.03142s^ 
a = -0.53733 - 0.60058s + 0.07863s^ + 0.13276s^ 
P = 0.75940 + 0.61356s - 0.43886s^ + 0.23942s^ 
7 = 0.10742s + 0.128005^ 



TV = 4.96269 + 1.54098s - 9.06376s^ + 4.94791s^ 
a = 0.05562 + 1.88660s - 2.943505^ + 1.04227s^ 
p = 2.79926 + 3.02991s - 4.148075^ + 1.914945^ 
7 = 0.85450s - 0.61016s2 



N = 4.25954 - 5.44309sc + 6.11031s^ - 3.13973s^ 



4 

a = -0.24144 - 1.07757sc + 1.52364s^ - 0.74308s^ 
P = 4.21265 + 0.25590sc + 0.98423s^ - 0.52893s^ 
7 = -0.04000sc + 0.08695s^ 



c 



N = 1.32443 - 1.41156S6 - 0.04809s^ + 0.79066s^ 
a = -0.88351 - 0.44818s6 - 0.60073s^ + 0.45526s^ 
P = 6.15221 + 0.46679sb - 0.50792s^ + 0.67006s^ 
7 = -0.00477S6 - 0.05503s^ 
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N = 0.23140 - 0.33644s + 0.162045^ - 0.02598s^ 
a = -1.36400 + 0.97182s - 0.02908.5^ - 0.431955^ 
(3 = 1.79761 + 1.57116s + 0.71847s^ - 0.68331s^ 
7 = 0.36906s + 2.39060s^ 

LO FFs for (p + p): 

N = 0.40211 - 0.21633s - 0.07045s^ + 0.07831s^ 
a = -0.85973 + 0.13987s - 0.824125^ + 0.431145^ 
f3 = 2.80160 + 0.78923s - 0.053445^ + 0.014605^ 
7 = 0.05198s - 0.04623s2 

N = 4.07885 - 2.97392s - 0.92973s^ + 1.23517s^ 
a = -0.09735 + 0.25834s - 1.522465^ + 0.77060s^ 
(3 = 4.99191 + 1.14379s - 0.85320s^ + 0.45607s^ 
7 = 0.07174s - 0.083215^ 

N = 0.11061 - 0.07726se + 0.05422s^ - 0.03364sj| 
a = -1.54340 - 0.20804s, + 0.29038s^ - 0.23662s^ 
p = 2.20681 + 0.62274sc + 0.29713s^ - 0.21861s^ 
7 = 0.00831sc + 0.00065s^ 

N = 40.0971 - 123.531S6 + 128.666s^ - 29.1808sg 
a = 0.74249 - 1.29639sb - 3.65003s^ + 3.05340s^ 
P = 12.3729 - 1.04932S6 + 0.34662s^ - 1.34412s^ 
7 = -0.04290S6 - 0.30359s^ 

Dp/p{x,fi^y. 

N = 0.73953 - 1.64519s + 1.01189s^ - 0.10175s^ 
a = -0.76986 - 3.58787s + 13.8025s^ - 13.8902s^ 
P = 7.69079 - 2.84470s - 0.36719s^ - 2.21825s^ 
7 = 1.26515s- 1.961175^ 
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NLO FFs for (p + p): 

N = 1.25946 - 1.17505s + 0.37550s^ - 0.01416s^ 
a = 0.07124 - 0.29533s - 0.245408^ + 0.165435^ 
p = 4.12795 + 0.98867s - 0.468468^ + 0.207505^ 
7 = 0.18957s -0.01116s^ (29) 



DP/P{x,fi^): 



Dfix,n') 



Dmx,ix^): 



N = 4.01135 + 8.67124s - 22.7888s^ + 11.4720s^ 
a = 0.17258 + 4.57608s - 9.64835s^ + 4.61792s^ 
(5 = 5.20766 + 7.25144s - 12.6313s^ + 6.07314s^ 
7 = 0.16931s -0.09541s^ (30) 

= 0.08250 - 0.04512sc - 0.00565s^ + 0.00900s^ 
a = -1.61290 - 0.38012sc - 0.06840s^ + 0.08888s^ 
(3 = 2.01255 + 0.63782se - 0.14146s^ + 0.06083s^ 
7 = -0.02958sc + 0.01130s^ (31) 

N = 24.2916 - 88.3524S6 + 93.1056s^ - 17.4089s^ 

a = 0.57939 - 0.80783s6 - 5.07200s^ - 2.45377s^ 
p = 12.1207 - 3.27370sb + 1.21188s^ - 5.50374sg 
7 = 0.14628S6 - 0.78634s^ (32) 



= 1.56255 - 1.48158s - 0.39439s^ + 0.51249s=^ 
a = 0.01567 - 2.16232s + 2.47127s^ - 0.93259s^ 
P = 3.57583 + 3.33958s - 3.05265s^ + 1.21042s^ 
7 = -0.84816s + 1.23583s^ (33) 
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Table 1: Xdf values and respective numbers of data points used for their evaluation for 
the various data sets shown in Figs. p]-|T2[ The selection criterion is a; > 0.1 in the case 
of e'^e~ annihilation and pt > f^o/^ otherwise. The cases where the experimental pt 
intervals contain values with p^ < fJ^o/C, are marked by asterisks. 



Type 



pp 



IP 



77 



Energy 
[GeV] 



133 
161 
172 
183 
189 



200 
500 
540 



630 

900 
1800 



165-251 



167-194 



10-30 
30-55 
55-125 
10-125 



Experiment 
DELPHI [H 



OPAL p3| 



DELPHI [22 



OPAL ^1 
DELPHI [H 

OPAL g3[ 
DELPHI [H 

OPAL g| 



OPAL p3| 



uAi m 



UAl [H] 

UAI y 

UA2 
UA2 0] 
UAl [H] 
CDF [g^ 
UAl [H] 
CDF [|i 



HI 3C 



Pt>2 GeV 
Pt>S GeV 
ZEUS 131 



OPAL 



Pt > 1.5 GeV 



OPAL |2 
Pt > 1.5 GeV 

OPAL PI 
Pt > 1.5 GeV 



OPAL §2\ 

Pt> 1.5 GeV 



1/2 



0.514 


6 


0.413 


11 


0.608 


6 


0.532 


11 


1.01 


6 


0.902 


11 


0.498 


6 


1.20 


11 


0.191 


11 


0.976 


14 


6.43 


14 


10.4 


14 



1.36 
7.68 
14.3 
27.0 
38.5 



12 
33 
2 
26 
19 



9.89 18 

59.9* 10 

14.8 10 

12.7 11 



2.58 
6.61* 
1.06 
97.8* 
2.13 
235* 
1.73 
70.3* 



1 



0.322 
0.317 
0.344 
0.353 
1.10 
1.04 
0.599 
1.64 
0.568 



6 
11 

6 
11 

6 
11 

6 
11 
11 



17.9 
5.13 
3.19 
2.30 
3.99 
1880 
2.04 
7.49 
6.56 



28 
28 
28 
6 
22 
40 
12 
40 
33 



0.467 26 

1.01 10 

0.564 10 

3.79 23 



5.38 
15.3 
1.70 
1.74 
2.33 
3.56 
3.25 
6.61 



12 

5 

14 
5 

14 
5 

14 
5 



0.498 
0.404 
0.360 
0.361 
1.26 
1.25 
0.877 
2.38 
1.57 



19.5 
24.6 
6.65 
12.4 
2.95 
6.89 
14.4 
26.1 



6 
11 

6 
11 

6 
11 

6 
11 
11 



48.4 35 

22.6 35 

18.9 35 

60.0 13 

21.0 29 

8360 43 

47.3 22 

27.2 47 

433 43 



4.15 26 

8.11 10 

4.11 10 

2.06 30 



15 
5 
17 

5 
17 

5 
17 

5 



19 



Table 2: Xbf values and respective numbers of data points used for their evaluation for the 
various pt distributions shown in Figs. |, and ||. The selection criterion is pt > 2/io. 



Type 



pp 



IP 



77 



Energy 
[GeV] 



200 
500 
540 

630 

900 
1800 



165-251 
167-194 



10-30 
30-55 
55-125 
10-125 



Experiment 



UAl 
UAl 
UAl 
UA2 
UAl 



m 



Hi 



CDF |9| 
UAl |25| 
CDF PI 



HI [|0| 

ZEUS m 



OPAL 



OPAL [0 
OPAL in 

OPAL m 



1/2 



0.976 
6.43 
10.4 
1.36 
7.68 
14.3 
27.0 
38.5 



9.89 
12.7 



2.58 
1.06 
2.13 
1.73 



4.93 
1.53 
1.44 
3.83 
354 
1.17 
4.08 
2.52 



0.416 
3.35 



2.18 
1.31 
0.110 
2.44 



9.68 
4.20 
3.32 
6.94 
756 
2.07 
9.41 
5.59 



2.65 
0.443 



2.80 
1.70 
0.467 
4.31 



No. of 
Points 



14 
14 
14 
12 
33 
2 
26 
19 



18 
11 



Table 3: Xdf values and respective numbers of data points used for their evaluation for 
the various x distributions shown in Figs. p!B|-P^. The selection criterion is x > 0.1. Data 



samples which did not enter the respective fits are marked by asterisks. 



Energy 


Flavour 


Experiment 


FF Set 


No. of 


[GeV] 






KKP 


13 


K 






BFGW ^ 




Points 


29 


uds 


TPC H 


0.178* 


0.159 


0.167* 


7 




c 




0.8 


76* 


0.911 


0.923* 


7 




b 




2.23* 


1.21 


1.14* 


7 


91.2 


all 


DELPHI 


1.28 


1.51* 


1.49 


12 






SLD i 


1.32 


0.370 


0.421 


21 




uds 


DELPHI [§ 


3.17* 


0.990* 


1.95 


13 






DELPHI [| 


0.201 


0.5 


88* 


1.00* 


12 




c 


DELPHI [| 


0.473* 


0.3 


88* 


0.401 


11 




b 


DELPHI [| 


2^ 


^.9* 


0.8 


87* 


1.03 


12 






DELPHI [| 


0.433 


9.14* 


8.74 


12 


189 


all 


OPAL m 


0.568* 


0.250* 


0.414* 


11 



20 




Figure 1: Normalized differential cross section {1 / a^otjdcy / dx (in pb) of inclusive charged- 
hadron production in e^e~ annihilation e^e'^ h'^+X as a function of scaled momentum 
X at CM energies y/s = 133, 161, 172, 183, and 189 GeV (from bottom to top in this 
order). The NLO predictions for ^ = 1/2 (dashed lines), 1 (solid lines), and 2 (dot-dashed 
lines) are compared with data from DELPHI |2^ (solid boxes) and OPAL ||2^ (open 
circles). Each set of curves is rescaled relative to the nearest upper one by a factor of 
1/10. 
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Figure 2: Differential cross section Ed^a/d^p (in mb/GeV^) of inclusive charged-hadron 
hadroproduction in pp collisions pp + X as a. function of transverse momentum pt 

at CM energies -y/i = 540 and 630 GeV, averaged over rapidity intervals \y\ < 2.5 and 
3, respectively. The NLO predictions for ^ = 1/2 (dashed lines), 1 (solid lines), and 2 
(dot-dashed lines) are compared with data from UAl The lower set of curves is 

rescaled by a factor of 1/10. 
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Figure 3: Differential cross section Ed^a/d^p (in mb/GeV^) of inclusive charged-hadron 
hadropro duct ion in pp collisions pp + X a.s a. function of transverse momentum pt 

at CM energies y/s = 200, 500, and 900 GeV, averaged over rapidity interval \y\ < 2.5. 
The NLO predictions for ,^ = 1/2 (dashed lines), 1 (solid lines), and 2 (dot-dashed lines) 
are compared with data from UAl [^. Each set of curves is rescaled relative to the 
nearest upper one by a factor of 1/10. 
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Figure 4: Differential cross section Ed^a/d^p (in mb/GeV^) of inclusive cliarged-liadron 
liadropro duct ion in pp collisions pp —>■ h^ + X as a function of transverse momentum pt at 
CM energy y/s = 540 GeV, averaged over rapidity intervals \y\ < 0.63 and 1 < \y\ < 1.8. 
The NLO predictions for ^ = 1/2 (dashed lines), 1 (solid lines), and 2 (dot-dashed lines) 
are compared with data from UA2 p7| , |28| . The lower set of curves is rescaled by a factor 
of 1/10. 
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Figure 5: Differential cross section Ed^a/d^p (in mb/GeV^) of inclusive charged-hadron 
hadropro duct ion in pp collisions pp + X as a. function of transverse momentum pt 

at CM energies ^/s = 630 and 1800 GeV, averaged over rapidity interval \y\ < 1. The 
NLO predictions for ^ = 1/2 (dashed lines), 1 (solid lines), and 2 (dot-dashed lines) are 
compared with data from CDF p9[. The lower set of curves is rescaled by a factor of 
1/10. 
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Figure 6: Differential cross section (Pa / dy dp\ (in nb/GeV^) of inclusive charged-hadron 
photopro duct ion in e^p collisions e^p —^h^ + Xa,sa, function of transverse momentum 
Pt at CM energy ^/s = 300 GeV, integrated over 7p-invariant-mass interval 165 < W < 
251 GeV and averaged over rapidity interval \y\ < 1 in the laboratory frame. The scattered 
positron is tagged, the maximum photon virtuality being Q^ax — 0-01 GeV^. The NLO 
predictions evaluated with the AFG @] photon PDFs for ^ = 1/2 (dashed lines), 1 (solid 
lines), and 2 (dot-dashed lines) are compared with data from HI 
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Figure 7: Differential cross section da/dy (in nb) of inclusive charged-hadron photo- 
production in e~^p collisions e'^p —>■ + X as a function of rapidity y at CM energy 
^/s = 300 GeV, integrated over 7p-invariant-mass interval 165 < W < 251 GeV and 
transverse-momentum intervals 2 < < 12 GeV and 3 < < 12 GeV. The scattered 
positron is tagged, the maximum photon virtuality being Q^ax — 0-01 GeV^. The NLO 
predictions evaluated with the AFG |0 photon PDFs for ^ = 1/2 (dashed lines), 1 (solid 
lines), and 2 (dot-dashed lines) are compared with data from HI The dot-dot-dashes 
and dotted lines correspond to the GRV p5| and GS pQ] photon PDFs. 
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Figure 8: Differential cross section (Pa / dy dp\ (in nb/GeV^) of inclusive charged-hadron 
photopro duct ion in e~p collisions e~p —^h^ + Xa,sa, function of transverse momentum 
Pt at CM energy ^/s = 296 GeV, integrated over 7p-invariant-mass interval 167 < W < 
194 GeV and averaged over rapidity interval —1.2 < y < 1.4 in the laboratory frame. The 
scattered electron is tagged, the maximum photon virtuality being Q^ax = 0-02 GeV^. 
The NLO predictions evaluated with the AFG ||3^ photon PDFs for ^ = 1/2 (dashed 
lines), 1 (solid lines), and 2 (dot-dashed lines) are compared with data from ZEUS 
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Figure 9: Differential cross section da/dpT (in pb/GeV) of inclusive cliarged-hadron 
pliotoproduction in 77 collisions e~^e~ e~^e~h^ + X as a function of transverse mo- 
mentum Pt at CM energy -^i = 166.5 GeV, integrated over 77-invariant-mass intervals 
10 <W <30 GeV, 30 < < 55 GeV, 55 < < 125 GeV, and 10 < W < 125 GeV 
(from bottom to top in this order) and rapidity interval \y\ < 1.5. The scattered electrons 
and positrons are antitagged, the maximum scattering angles being 6' = 33 mrad. The 
NLO predictions evaluated with the AFG |Q photon PDFs for ^ = 1/2 (dashed lines), 
1 (solid lines), and 2 (dot-dashed lines) are compared with data from OPAL [^. The 
lowest two sets of curves are rescaled relative to the nearest upper one by a factor of 1/10. 
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Figure 10: Differential cross section da/d\y\ (in pb) of inclusive charged-hadron photo- 
production in 77 collisions e~^e~ — > e~^e~h'^ + X as a function of rapidity \y\ at CM 
energy ^/s = 166.5 GeV, integrated over 77-invariant-mass intervals 10 < < 30 GeV, 
30 < < 55 GeV, 55 < W < 125 GeV, and 10 < W < 125 GeV (from bottom to top 
in this order) and transverse- momentum interval 1.5 < pr < 15 GeV. The scattered elec- 
trons and positrons are antitagged, the maximum scattering angles being 6' = 33 mrad. 
The NLO predictions evaluated with the AFG ||3^ photon PDFs for ^ = 1/2 (dashed 
lines), 1 (solid lines), and 2 (dot-dashed lines) are compared with data from OPAL ||3^ . 
The lowest two sets of curves are rescaled relative to the nearest upper one by a factor of 
1/10. 
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Figure 11: Differential cross section da/dpT (in pb/GeV) of inclusive cliarged-hadron 
photoproduction in 77 collisions e~^e~ — > e~^e~h^ + X as a function of transverse mo- 
mentum Pt at CM energy -y/s = 166.5 GeV, integrated over 77-invariant-mass interval 
10 < W < 125 GeV and rapidity interval \y\ < 1.5. The scattered electrons and positrons 
are antitagged, the maximum scattering angles being 6' = 33 mrad. The NLO predic- 
tions evaluated with the AFG |0 (solid line), GRV [^] (dot-dot-dashed line), and GS 
Il36| (dotted line) photon PDFs for ^ = 1 are compared with data from OPAL [^. For 
comparison, also the direct (dashed line), single-resolved (dot-dashed line), and double- 
resolved (dash-dash-dotted hne) components of the AFG prediction are shown. 
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Figure 12: Differential cross section da/d\y\ (in pb) of inclusive charged-hadron photopro- 
duction in 77 collisions e+e~ e'^e~h^ + X as a function of rapidity \y\ at CM energy 
\/s = 166.5 GeV, integrated over 77-invariant-mass interval 10 < ly < 125 GeV and 
transverse- momentum interval 1.5 < < 15 GeV. The scattered electrons and positrons 
are antitagged, the maximum scattering angles being 6' = 33 mrad. The NLO predic- 
tions evaluated with the AFG |0 (solid line), GRV [^] (dot-dot-dashed line), and GS 
||36| (dotted line) photon PDFs for ^ = 1 are compared with data from OPAL [^. For 
comparison, also the direct (dashed line), single-resolved (dot-dashed line), and double- 
resolved (dash-dash-dotted hne) components of the AFG prediction are shown. 
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Figure 13: Normalized differential cross section {1 / a tot) d<^/dx (in pb) of inclusive charged- 
hadron production from b, c, light, and all quarks (from bottom to top in this order) in 
e~^e~ annihilation e'^e" ^ /i^ + X as a function of scaled momentum x at CM energy 
^/s = 29 GeV. The NLO predictions for ^ = 1 based on the KKP g (sohd lines), K 
||T9| (dotted lines), and BFGW ||2^ (dashed lines) FFs are compared with data from TPC 
p5|| . Each set of curves is rescaled relative to the nearest upper one by a factor of 1/5. 
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Figure 14: Normalized differential cross section {1 / aiot)d(y / dx (in pb) of inclusive charged- 
hadron production from 6, c, light, and all quarks (from bottom to top in this order) in 
e^e~ annihilation e^e^ ^ /i^ + X as a function of scaled momentum x at CM energy 
= 91.2 GeV. The NLO predictions for ^ = 1 based on the KKP H (sohd lines), 
K [|T9| (dotted lines), and BFGW (dashed lines) FFs are compared with data from 
DELPHI (yields of 1991-1993 (triangles) [| and 1994 (circles) 0) and SLD (squares) §. 
Each set of curves is rescaled relative to the nearest upper one by a factor of 1/5. 



34 




Figure 15: Normalized differential cross section {1 / a tot) d<^/dx (in pb) of inclusive charged- 
hadron production from b, c, light, and all quarks (from bottom to top in this order) in 
e~^e~ annihilation e'^e" ^ /i^ + X as a function of scaled momentum x at CM energy 
^/s = 189 GeV. The NLO predictions for ^ = 1 based on the KKP (solid hues), K [0 
(dotted lines), and BFGW ||20[ (dashed lines) FFs are compared with data from OPAL 
[]23|| . Each set of curves is rescaled relative to the nearest upper one by a factor of 1/5. 
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